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I. INTRODUCTION 


The shielded suspended-substrate line is a transmission medium useful for 
radar and microwave circuits in the Ka band frequency range, 30 — 40 GHz. 
In order to utilize this transmission medium in the construction of microwave 
circuits and filters, it is necessary to have valid circuit models for typical 
discontinuities such as the series gap in line, open-ended stub, and a 
discontinuous change in width. These discontinuity characteristics may be 
deduced on the basis of a calculation of their scattering coefficients for 
propagating waves in the medium. The approach chosen in eet 
however, is that of placing the selected discontinuity in an open-ended strip 
resonator in the transmission medium. An appropriate circuit model is 
adopted for the given discontinuity, and its circuit elements are then deduced 
on the basis of the perturbation of the resonance frequency of the strip 

oS ee 

The boundary-value problem associated with the microstrip resonator 
structure has been approached in a rigorous manner based on a full-wave 


rman ae ee 


analysis which utilizes a process of solving the electromagnetic (EM) boundary 
value problem with inclusion of all the field components. “An analysis : V4 2 


presented by Itoh and Uwano was employed in the present work (Refs. 1, 3]. 





The characteristic equation for resonant frequencies was obtained by use of 
Galerkin's technique applied in the spectral or Fourier transform domain. The 
set. of algebraic equations among the transformed quantities thus produced is 


called the Green's function relations in the transform domain. 








Il. THEORETICAL ANALYSIS 


A. FIELD EQUATIONS 

The microstrip resonator to be analyzed is shown in Fig. 1. A 
rectangular strip of width 2w and length 2¢ is placed on the suspended 
substrate. The sides and the top of the structur’ are surrounded by 
conducting shield walls. Thus the entire structure is considered to be a 
suspended-strip resonator located in a partially filled waveguide with end 


walls. 





Figure 1. End and top view of Suspended-Substrate Microstrip Line 


It is assumed that the thickness of the conducting strip is negligible and 
that all the media and conductors are lossless. For simplicity, the strip is 
assumed to be symmetrically located, although the present method of analysis 


can be easily extended to more general cases. 














The fields existing in the structure shown in Fig. 1 are the superposition 
of TE (to z) and TM (to z) fields, which can be expressed in terms of two 
types of scalar potentials $¢(x,y,z) and $4(x,y,z), where the superscripts e and 
h denote electric and magnetic, respectively [Ref. 1]. 























a°,° 
E,;(xy,2) = k;%,° + — po (1a) 
3o;" 
H (xyz) = k,7o;5 + (1b) 
022 
2d .e h 
E,;(xy.2) = ae — jp; “ (1c) 
Hys(4y.2) = ae je, a (14) 
where i = 1, 2, 3 designates the substrate or the air region, as shown in 
Fig. 1. 
k,? = ka? = kg? = ppegu” (1e) 
ky? = pyégu” = pe ko? (1f) 
& = €3 = & (1g) 
bo = €€9 (1h) 
th = Hs = Uo (14) 
Hy = Mello = Mor (Hy = 1) (1)) 


where w is the operating frequency and €) and pp are the free-space 


permittivity and permeability. 





RR Pee noe ee. eee a 





B. BOUNDARY CONDITIONS 
Now applying boundary conditions at the bottom in region 1, tangential 


electric field and normal magnetic field components must be zero. 


7] 

E,,(x,0,z) = 0, —H,,(x,0,z) = 0 (2a) 
dy 
a 

E,,(x,0,z) = 0, pps) = 0 (2b) 


At the interface between regions 1 and 2, the tangential field components must 


be continuous. 


At y = h;: 
E.1(x,h,,2) = E,,(x,h,,2) (3a) 
Ex4(X)by,2) = Eyo(x,h,,2) (3b) 
H,,(Q%h,,2) = H9(x,hy,2) (3c) 
H,,(x,hy,2z) = H,,(x,h,,z) (3d) 


At the interface between regions 2 and 3, tangential electric field components 


must be continuous. 





At y = hy: 
E,o(xjho,2) = E,3(X;h9,2) (4a) ; 
Byo(x,hy,2) = E,,3(x,hy,2) (4b) | 
Also, the electric fields at y = hy will exist only on the air/dielectric 


interface, and can be expressed as, i 





E,(x)e”, a> |x| 2 w, : 

E,9(x;ho,2) = ; (4c) j 
0, otherwise 

E,(xe™, a> |x| 2 w, 

Ex3(Xsho,2) = : (4d) 
0, otherwise | 








Similarly, tangential magnetic fields at h, must be discontinvous by the 


corresponding surface current densities. 


Je”, Ix] <w, 
H9(x,ho,2) - H(x,h),2) = . (4e) 
0, otherwise 
J,(xje™, |x| < w, 
Hyo(x,ho,z) — Hyg(x,hy,z) = ; A (4f) 
0, otherwise 


At the top in region 3. tangential E-field components must be zero. 








E,3(Xsh3,z) = 0, (x,hg,z) = 0 (5a) 
dy 
OH, 

E,.3(Xshg,z) = 0, ‘ (x,hg,z) = 0 (5b) 
y 


The final boundary conditions occur at x = + = where the tangential E-field 


components must be zero in all regions. 


At x = ao 
E,i(tsoy,2) = 0 (6a) 


Ey(@-Soy2) = 0 (6b) 


aI] 








etek ie te 





C. SPECTRAL DOMAIN APPROACH. 

The boundary value problem is solved in the Fourier transform or spectral 
domain. This leads to the Green's function equations associated with the 
structure. This analysis was first carried out by Itoh and Mittra and was 
successfully applied to a number of problems in microwave integrated circuit - 
structures [Ref. 7]. 

1. Potential Equations 

The quantities in Eqs. (1) are transformed into spectral domain via 
the Fourier transform; 

6¢8(n,y,6) = i dz ie dxh?"(x,y, zelkn Xejbz (7a) 


n= 1, 2, 
where # is the Fourier transform variable, k, is the discrete transforru variable 


defined by, 











(n-—-)r 
k, = ———— for E, even, —H, odd(in x) modes. (7b) 
a 
nen 
k, = for E, odd, —H, even(in x) modes. (7c) 
a 
The transforms of field quantities are now, 
E.;( (n,y,8) = k, 3 2 6;° (8a) 
H,,(n ny,8) = kd," | (8b) 
it 
B(ny,B) = ky Ab? - joy; (8c) 
_e& 
H,;(n,y,8) = -k 4," na jwe; (8d) 
kg = kp + M, (P= -j8, k? = wy,e)) (Se) 
kg = ke - P= pteegQ? - FP si 


where i = 1, 2, 3, for each of the three regions defined in Fig. 1. 














The transforms of scalar potentials satisfy, 


d?perh 
: x 92g = 0 (9a) 





1? = 737 = ky? — kei 3 =k,’ + #? - 
k;? a kj = k,? + p? = € Uk? 


The solutions for these homogeneous differential equations are well known 
and can, be described in the general form as below, 


6; = C,coshyy + C,sinhyy. (9b) 
i 1 2 


Now the boundary conditions (2), and (5) are applied. The cosh term 
vanishes in electric field and the sinh term vanishes in magnetic field in 
regions 1 and 3. The potentials are: (see details in Appendix A), 

Region (1) 


g,° = Asinhy,y (10a) 
g," = Beoshyy (10b) 
Region (2) 
$,° = Csinhyy + Deoshyy (10c) 
g," = Esinhyy + Feoshyy (10d) 
Region (3) 
$,° = Gsinhy,(y — hg) (10e) 
$5" = Heosh»,(y — hg) (10f) 
where 
ye = ky? fi ala (10g) 
o" even 
onde a - (10h) 














2. Green's Function 
In conventional space domain analysis, the structure may be analyzed 
by first formulating the following coupled homogeneous integral equations. 


The equations must then be solved for the unknown propagation constant £. 


f [Z2(x-x' y)J,(x') + Z,.(x—x',y)J,(x')]dx' 


f (Z,..(x-x',y)J,(x') + Z,.(x-—xy)J,(x')]dx! 


E,(x) (11a) 


E, (x) (11) 
where E, and E, are unknown electric fields on the boundary at y = hy, J, 
and J, are current components on the strip at y = h,, and the Green's 
functions 2; (i,j = z,x) are functions of £. 

The integration is over the strip where E,(x) and E,(x) are zero, as the 
strip is assumed to be perfectly conducting. In the spectral domain 
formulation, however, the following algebraic equations, instead of the coupled 
integral equations, are obtained. These equations are the Fourier transform 


equivalent of the coupled integral equations. 
Zo(Kysho)J,(k,) + Z.x(Kysby)J,(k,) = E,(kpshy) (12a) 
Zy2(Kysho)I (Kp) + Zyx(Kysto)J, (ky) = B,.(Kysh) (12b) 


where quantities with tildes (~) are Fourier transforms of corresponding 
quantities. The Fourier transform is defined as in Eq. (7). 

The right-hand side of equations (12) is not zero everywhere because the 
Fourier transform requires integration over all x, not only over the strip. The 


equations contain four unknowns J,, J,, E 


» E, and E,, with unknown £. 
However, E, and E, will be eliminated in the solution process based on the 


Galerkin procedure. 
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In the spectral domain, the boundary conditions (3), (4) are now given by 


. the following equations. 


at y = hy 
BE. = Ey 
Ey <= Eo 
Hy = Hy» 
Ai, a Hy 
at y = hy, 
E. = Ej 
Exo = E,s 
Hyp — Hyg = J, 
Hy = Hy a J, 


where J, and J, are Fourier transforms of unknown current components J,(x) 


and J,(x) on the strip at y = hy. 


1] 








aR canteen 








When the fields are expressed in terms of the potentials, these can be put 
into the form: [Ref. 3] 








i - Hig 2 i (13) 
FE, Lx. Zxx x hd 
where 
, ee casa 
Lag = oP (PZ, + k,°Z,] 
Laz = 2, = - ——{, - 4] 
k + 
be = - ik", + 2] 
k 2+" 
‘, = Yy2Ct, + Yy1Ct, 
Ct,Ct. + CtyCty Re + CtyC tg ft + a 
i, = %a0t, + Yat 
Var VqQCtCty + YerIegCtiCts + Y22%gCteCts + 7,3 ; 
For 7, 2 0 


Ct, = coth7,d, Ct, = coth yt, Ct; = coth 7h, 
For 7, < 0 





Ct, = -jcoty,d, Ct, = -jcotyt, Cts = -jcotygh, 4 


i =Vk 2 + P- lyeiur 


%% a 
Vi = 7 Ng oe . 
jue; jun; 








The quantities Z,,, Z,,, Z,, and Z,, are actually the Fourier transforms of 


dyadic Green's function components. ei 








3. Characteristic Equations in Matrix Form 
Two equations in (13) contain four unknowns E,, E,, J, and J, 
However, the first two unknowns E, and E,, can be eliminated by applying 
Galerkin's method in the spectral domain. The first step is to expand the 
unknown J, and J, in terms of assumed basis functions J,,, and J, with 


unknown coefficients c,, and d,,. 


J,(n,8) = Dog) am(2, ) (14a) 
5,(0,8) = B.4y5_q(0,8) (140) 


The basis functions J,,, and J,,, must be chosen to be the Fourier transforms 
of space-domain functions J,.(x,z) and J,,,(x,z) which are physically realistic, 
and which are zero except for the region |x| < w and |z| < @ Now, 
substituting (14) into (13) and taking inner products of the resulting equations 


with the basis function J,, and J,, for different values of i yields the matrix 


equation, 
Bier, No. Sip ie, VG ot 
f Jeidze endo + JriZex % In) xml = 0, (15a) 
ee Oe ,M 
J i2ee Beg Jn + IgiZex > bd manld3 = 0, (15b) 
es ee ,N 


The right-hand sides of (15) are zero by virtue of Parseval's theorem, because 
the currents J,;(x), J,;(x) and the field components E,(x,hy), E,(x,h)) vanish 
in complementary regions of x. For example, if the inner product of E,, on 


the left-hand side of (13) with J,;(k,) is taken, one obtains 


fi Ji(k,)E,(k,)dk, = anf J,i(x)E,;(x)dx = 0 


In the above, J,;(x) is zero outside the strip, and B,(x) is zero on the strip. 
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Therefore, the final boundary condition is now used. 


Equations (15) will be expressed in matrix form as follows, 


N M 
‘Kr ( 11) “1 (192) as 
PAST Cn + eK io dy 0, Gs i 2, sea N. (16a) 
x 291 . 252 
Kim Cy + RST y dn = 0, : (16b) 


It 


i= 1,2, -,M. 


where from the definition of the inner products associated with the Fourier 


transform defined by (7), the matrix elements are 


151 ft % 7 TF 

K{i) (ky) = 3 f J,j(0,8)Z,Jom(0,)d8 (17a) 
o oO eae 

Kia (Ko) = BJ, xilsP)ZerJxm(s8)4 (17b) 
wo oO, sad 

Kin? (ko) = BJ Ixi(tB)ZxaSom(ns8)d8 (17) 
oO 0, Petes fd 

K (29?) (kp) = Py 0 J §(1s8)Z 55 xq (0,8)d8 (17d) 


A homogeneous system of equations is thus obtained in terms of the 
unknown coefficients c, and d,. In order that c,, and d, have nontrivial 
solutions, the determinant of the matrix must be zero, and hence the resonant 
frequency w is determined for the resonator and discontinuity assumed. 

Equations (16) are now solved for the wave number ky by ‘setting the 
determinant of the coefficient matrix equal to zero and by seeking the root of 
the resulting characteristic equation. The resonance frequency of the 
suspended-stripline resonator is derived from the obtained value of kp. 

The accuracy of the solution can be systematically improved by increasing 
the number of basis functions (M+N) and by solving larger size matrix 


equations. However, if the first few basis functions are chosen so as to 





| 
| 
| 








approximate the actual unknown current distribution reasonably well, the 
necessary size of the matrix can be held small for a given accuracy of the 
solution, resulting in numerical efficiency. Hence the choice of basis functions 


is important from the numerical point of view. 
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Jnu(z) = b (20a) 
1 


ae (f+), ses e4=m-s 
J,o(z) = . 
“37, sin - 5, gsz6p.,4=7,-8 


With a gap present assume e = 0, and J, = CJ, + CoJ,9, with distributions 


(20b) 


as shown in Fig. 3. The transformed distributions are: 
Jay = Jny(B)-Joy(ky) (21a) 
Joa = Jza(B)*Iex(kp) (21b) 
J.x(k,) is same as in Fig. 2. 
 lPPs = @-IPP2 ~ visin( 6g) 
Alp, + Po — 2g) 


Jzo(8) = 
: Re 
15 at Bg aa 2) a OS 
—le 2 +e 
4 a+ Bd, 7 BPo 
38t ithe, et As «st alr-Ap Az - As 
J & J = 
ef e été a2 | (22) 
t + Bd, T - Bd, 
From Eq. (15a) and with i = 1, 
JaiZerm Jom = Jng(1,B)ZonC 54 + Ja(0,8)ZyeCod a0 (23a) 
and, with i = 2, this takes the form 
sy An ae Joo(0,8)Zp_Cy5 94 + Jz9(0,8)Z_qCo) 99 (23b) 


These equations can be represented in the matrix form: 











K1l K12 | K 
= (24) 
K21 K22 | le, 0 
where, 
K1l = Sf 5,4(n,8)2Z,,5,1(n,6)46 K12 = Sf F,4(0,8)Z,,I29(n,8)4B 
K21 = Sf F.o(n,)Z,.5,1(n,8)4B K22 = Df 5,9(n,B)Zood29(n,8)d9 
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B. BOUNDARY VALUES 

The actual computation was carried out using the shield and line 
dimensions shown in column i, in Table 1, which is typical data from a 
suspended stripline filter problem initiated by the NOSC (Naval Ocean 
Systems Center) Microwave Laboratory in San Diego. This configuration is 
used’ in the frequency range of 30 - 40 GHz. 

The data in columns 2 and. 3 were used in analysis of the microstrip line 
with a gap, to compare with published data {Ref 5]. 

Data 2 uses frequency range of 85 — 102 GHz 

Data 3 uses frequency range of 31 — 40 GHz 

TABLE 1. CONFIGURATION DATA 


(dimensions : mm) 





Symbols for Fig. 1 
.a = Shield width | 
d = Height of lower air layer 
t = Height of dielectric substrate layer 
h = Height of upper air layer 
w = Microstrip line width 
g = Gap in microstrip line(center) 
= Er = dielectric constant of substrate layer. 
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4, Analysis of Resonator with Gap Present 
a. Resonance stub lengths at given frequency pairs 
The moment method calculation (Appendix D. part A) was 
carried out with various assumed stub lengths. The results are shown in 


Table 2 and Table 3. 


TABLE 2. RESONATOR STUB LENGTHS WITH DATA 2 (SEC. II-B) 


Units 
length : mm 
frequency : GHz 


oe 


0.05 1.211647 1.075264 0.959969 
(1.230606) (1.090897) (0.973572) 


1.259331 1.121864 1.005066 
(1.278259) (1.137848) (1.014435) 


1.285843 1.147675 1.030686 
(1.304749) (1.163574) (1.044349) 
1.300917 1.163193 1.045738 
_ (1.319977) (1.179146) (1.059464) 
0.25 1.310349 1.172264 1.055008 
(1.329392) (1.188179) (1.068817) 
1.316086 1.178024 1,060890 
(1.334999) (1.193985) (1.074570) 
































c. Effective dielectric constant 
The calculation of the effective dielectric constant used the 6 values 


which were computed above (Table 4). 


d. Characteristic impedance Z, [Ref. 6] 
The computation of characteristic impedance uses an empirical 
calculation for the quasistatic valuc, Z)(0) (Ref. 5] which has been confirmed 
to be accurate [Ref. 8]. This value is then corrected to the wanted frequency 


f by use of €,¢(f) in the expression: 


Zo(0) 
Zo(f) = wea (Appendix D. part B. data file) 
Vv eff 
(1) for0<w< =; 
Z, = 60[V . Rlnfe2 a ee (2) 
Z 
tyes 
Vv Cerf 


V = -1.7866 - 0.2035(--) + 0.4750(-+) 
R = 1.0835 + 0.1007(++) — 0.09457(-+) 





(2) when 2 <w<a; 





2 
Z, = 120z[V + R{--— + 1.3930 + 0.6670-én(~ + 1.444)}] ) 
Z, } 
Zo = ——_. | 
V = -0.6301 ~ 0.07082(L) + 0.24704) 
R = 1.9492 + 0.1553(-7.) - 0.5123(-4) “ 
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IV. DERIVATION OF THE PI-CIRCUIT 
REPRESENTATION OF THE DISCONTINUITY 


A. EQUIVALENT NETWORK 

The gap discontinuity is represented by a capacitive Pl—network. This 
circuit is terminated on each side by an open-ended stub. The PlI-circuit 
components are evaluated in terms of the admittance matrix, [Y;,], of the 


two-port discontinuity, as shown below. 





Figure 11. PI—Network 


The admittance of an open line segment is: 


Yin = Y AB = JYotanpt, (€ = strip line length) 





From the 2-port admittance matrix, 
T,-= YuVy + YyVo 
Ip = YorV, + YooVo 
I, Vo 
re th Nig 
1 1 - 


Deere er ses 





36 








| — = Yo,— + Yx 
| Vo Ne 
aa ae ee Yo9} 
Vo V2 You 
: Vo You 
V; I, 
Fe 
Vo 
— = Yu t Yy 
1 2. 
fe os 
But at resonance 
V, 1 Vv, 1 
ne ee |e 2 
YpYor 
-Y, ~ Yu = 
: ~Y¥y- ¥en 
[Yu +, } [Yn é 4) -~Y2=0 (25) 


This equation is used to calculate the equivalent circuit capacitances. The 





admittance—-matrix values for the equivalent circuit parameters are as shown in 


Fig. 5. 





4 








a Figure 12. Admittance network 
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B. SUSPENDED STRIPLINE GAP CAPACITANCE 

Figure 13 is the equivalent network of Fig. 1 which shows the 
configuration of the gap discontinuity in a suspended stripline. This network 
is enclosed in the central box in Fig. 14. In the computation procedure, stub 
lengths ¢, are assumed and the corresponding frequency of resonance is found. 
Then a slightly different length 4 is employed to find the new resonance 
frequency. The resulting two solution values provide the necessary data for 


the calculation of the discontinuity capacitances, C,,, and C,,, 


Yui + Yi2 = jwCpar 





Figure 13. Network parameters 


at eo 


= -jYotanfl, = -jYotan£l, 


se 





Figure 14, Analysis procedure 


(Yu + ¥ Yn + Y,) = Yp (26a) 
t ' 
(Yu + Y,(¥n + Y,) = YS (26b) 
9 t Pe \ - 
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' of 
—  YYe-¥4Y%, YC) + ¥,(6) 
ig BS 
Y, +Y. -¥Y, =¥ 2 


TABLE 5. GAP CAPACITANCES WITH DATA 2 (SEC. a 














2 
io 


TABLE 6. PARALLEL CAPACITANCES WITH DATA 2 (SEC. III-B) 


cev[_* [ «| @ 
0.25 14.0813 14.6777 
0.3 13.7098 14.2771 


* In Tables 5 to 8, frequencies f are given in gigahertz, and the 
tabulated capacitances are in femtofarads. 
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TABLE 7. GAP CAPACITANCES WITH DATA 3 (SEC. III-B)” 


0.7176 
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C. RESULTS 

Tables 5 to 8 summarize the series and parallel equivalent-circuit 
capacitances of the gap discontinuity over the range of calculated dimensions 
and frequencies. This data is also shown in graphical form in Figs. 15 to 18. 


As may be seen in these Figures, both the series gap capacitance (C and 


gap) 
the shunt parasitic capacitance (C,,,) were found to decrese monotonically 
with gap length, at all frequencies. The plots of capacitance vs. frequency 
show, however, that while the parallel capacitance rises with frquency, for all 
system dimensions computed, the series gap capacitance Ceap rises to a 
maximum near the center of the K-band region, but C,,, shows a trend 


toward minimum values at mid-w-band frequencies. 
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Figure 16. Capacitances at W—-band 
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V. CONCLUSIONS 


This thesis follows the methods originated by Itoh for the analysis of the 
suspended stripline resonator (Refs. 1 and 3], and for the equivalent circuit of 
the gap discontinuity (Ref. 4]. Some changes have been introduced in these 
approaches, however. A fully-shielded enclosure was employed, allowing the 
use of a finite Fourier transform in two coordinate directions. This change 
reduced the computation time significantly, while reproducing the results of the 
integral transform along the line axis. 

The perturbed-resonator technique permitted the use of the strip current— 
density distributions suggested by Itoh [Ref. 1]. These are known to give 
accurate results in the Galerkin method. A calculation was made employing 
the resonator configuration used in [Ref. 5], which is the only known literature 
reference on the suspended stripline gap. 

Table 5, Table 6, Fig. 16 and Fig. 18 shows the resulting parameters at 
W-band. Table 7, Table 8, Fig. 15 and Fig. 17 show the parameters at Ka- 
band. It is evident that all these capacitances are strongly frequency- 
dependent. Comparing these data with data from (Ref. 5], we see that the 
gap capacitances found in the present work range from 0.1 - 1.5 femtofarads. 
These gap capacitances are smaller by an order of magnitude than those found 
by Rong and Li (Ref. 5], but the accuracy of the latter work has not been 
confirmed. When the gap is widened, the capacitance changes at a rate 
comparable with that in (Ref. 5]. The parallel capacitances found here range 


from 13 — 70 femtofarads. This parallel capacitance shows an inverse 
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dependence on gap length, unlike that of (Ref. 5]. The reason for this 
discrepancy is not known, although it is felt that the capacitance variation 
found in the israel work is correct. 

The two-frequency calculation for each capacitance used here introduces 
some error, but this should be small for small computed frequency differences. 
In any event, the circuit-equivalent capacitances found in this way should be 
accurate at a frequency intermediate between the two frequencies of the pair 
employed in each case. 

Recommendation: In attempting to use Itoh's method of perturbed 
resonators having the same frequency, the determinant of Eqs. (17) did not go 
to zero. The reason for this failure should be determined, so that an ideally 
accurate computation of equivalent-circuit parameters at a well-defined single 


frequency can be made. 
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APPENDIX B 
FOURIER TRANSFORMS OF THE ASSUMED CURRENT DISTRIBUTIONS 


The evaluation of the Fourier-transformed forms of the strip current 


densities is shown in detail below: 


A. CONTINUOUS STRIP 


The coordinate form of the current distributions in Fig. 2. 











1 
J,(x) = —1 + |-*- il -w<xiw (B-la) 
2w 
1 
J,(z) = -s cosr-, <2 (B-1b) 
Jg(x) = —— sin as ; -w<xiw (B-ic) 
w 
Z : 
J,(z) = ve ~<2¢ (B-1d) 


The Fourier transform is defined to be: 


Ha) = f oxelMax (B-2) 





The Fourier transforms are calculated as shown: 


i(n) = face ax, (B-a) | 
5,(6) = fseyele, (B-3b) | 
jn) = f sje as, (B-3c) ; 
1A) = f S4ceay, (B-3d) 











t 


The Fourier transforms of the basic current distributions are: 


1,8) = f 4- cos elas, 46260, 
~ @ 7 
J,(n) = f sin X_elknXgy, -w <x <w, 
0 
0 ‘ 
5,(8) = f —a2-22ae, 6266 
(8 = f are 


The Fourier transforms are calculated as shown below: 


From Eq. (B-4a) 


(B-4a) 
(B-4b) 
(B-4c) 


(B-4d) 


x © ‘ w ‘ 
j(n) = fel rn Ea alia = 1. f E ‘ ray aii 





wW 1 ; a 
Ww 
= ae [f eketex + Sf Ixl°elo*as] 
mM ie =a 
The first term of Eq. (B-5a) is evaluated as follow: 
1 w 
f ebkiXay = __alk,x 
“W jk, X=-w 
= —elkn¥ _ etn) 
jk, 
1 
= ie (cosk,w + jsink,w —- cosk,w + jsink,w) 
n 
pane SV 
= Bau 


The second term of Eq. (B-5a) is evaluated as follow: 


Ww k : 
f Ix[8e*o%ax 
-W 


Use fudv = uv — fvdu 
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(B-5a) 


(B-5b) 


(B-5c) 

















For positive in Eq. (B-5c) we have: 












































f "delkaX ay 
0 : 
‘ jk,x a jk,x 
= x3 = f 3x2 — dx 
jk, 0 jk, 
x3¢ Jkpx 3 , x2edknx w edKnx 
= — —|—__ - f 2x——<x| 
jk, jk, jk, 0 jk, 
x3e Jkax 3 , x2edknx 2 - xedknx w edknX 
Saini Lalas ately pa 
J ky jk, J ky jk, jk, 0 jk, 
x3e Jknx 3x2e En 6 , xedKnX oJ knx W 
2 zs . ; eee }| 
jk, kn k jk, k x=0 
we jk,w Swie jk,w 6 wedknW jk,w F 
ee ee St 
ik, ky ky iky ky ky ; 
wieJknW —gwedknW ye dkaW gg dKnW 6 
SS (B-5d) 
k, Ke ki k,4 Kas 
For negative x in Eq. (B-5c) we have: 
- f ” <del gy 
-w 
x8elknX genx gg xe@dknX = edn |g 
A ee 
jk, ky? k,” jk, k,* xX=-w 
ey ae eee es a | , 
Kee jk, k, Ki jk, ky” : 
6 wie JKyW swe Jka gwe ska ged knw 
= + jp 4 — -j - — _ ) 
ky kn, k, ky ky 











Now we add. up together Eq. (B-5b) Eq. (B-5d) and Eq. (B-5e) into Eq. (B-5a) 





























; 1, 2sinkw 1 weelkn™ wed Kn 
J,(n) = = + —.| -j————“—- = ee 
2w k, w k, ky 
ewelknW — gelKn® 4g 6 wee Jk gw2e7S Kaw 
urge ee Se ee ek ee = ee 
ky k, ki ky ky ky 
ewe Jaw ge Skaw 
eS 
k,? ka! 
1 > 2sink,w 1 waes knw wre J knw swed kaw 
= _— zt —| ee ee eee 
aw t k, w ky k, ky? 
sweedKn™ wed tn® — gwekn™ gedkn™ gedkn¥ ap 
+ ——_—__ + j ____- j___- aeons +—}| (B-5f) 
ky” k,° k,° kn Ky’ ky4 ; 
Now we consider Euler's equation 
ekaW _ eIknW = oisink w 
kaw 4 eIn¥ = 208k, w 
. 1 > 2sink,w 1 ,  w(2jsink,w) 3w?(2cosk ,W) 
\n = an eee = |43——— + ar rrr er 
2w ki, w k, ki 
_ 6w(2jsink,w) 6(2cosk,w) 12 
ee a ay 
ky ky ky 
sink , W 1 2w°sink,W 6w’cosk,w  12wsink,w 
ey Cy (er + oe 
k,W 2w4 ky ks his 
12cosk,w 12 
k,4 k,! 
_ 2sink \w 3cosk ,w 6sink,w 6cosk ,W 6 
k,W k,?w? k, ow? k, tw! k, {wt 
dsink,w 3 2sink ,w 2(1 - cosk,w) : 
= + ae {cosk yu of | -5g) 
k,W k,-w" k,W k,°w" 




















= - | —dz 
2 je -£ jp 
2 ip mF zat 


20 ' jp (ig? ke 
= {ell 5 rity 4 ell erist) 
2° ‘ jp 
t 1 
= salty : sina} 
ae ical 7 eiaedl (3-8) 
pe (Be)? 


Finally the current distributions are found to be: 
From Eq. (B-5f) and Eq. (B-6) 
J,(kys8) = J,(n)-5,(8) 
2sink ,w 3 
-—-+ 
k,W k,?w? kW kw? 
From Eq. (B-7) and Eq. (B-8) 
J<(ky18) = Ja(n) +5 4(8) 
2msink , W ( [ cos Be sinGe } 
wm —(k,w)? Be (Be)? 
2rsink , W cos fl sa 


Qsink,w  2(1-cosk,w))7 4acosfe 
rae ae oT 
7?-4(0) 








{cosk,w . 








(B-9) 





= 


(B-10) 
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B. RESONATOR WITH GAP 
Asume J, = 0 
J, = C31 + Codze 
Jor = Jny(z)-J,,(x) 
Jaa = Izo(Z)*Izx(X) 
In both open-ended line segments of same length, current J,,(z) is chosen 


to have the form shown in Fig. 19. 





Figure 19. Current distribution of J, 
J(2) = b, |x| $ a, |2| cae A 


Normalization of J,,(z) leads to: 


(p, — 8) + (pp - g)b = 1 
1 
———————— 

Pi + Po — 2g 
1 
J,,(2) = b = —————— (B-11a) 
Py + Po - 2g 


Using J,o(z) from Fig. 3, we have: 
Let d, = p, - : a = Pp - 
- 577 sin( 4 Hy -Pp <2 <-g 


a sin( i a, ge<z<p 


Jzo(2) = (B-11b) 


The Fourier transform from Eq. (B-1la) is: 


=f J,,(2)0 az 
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APPENDIX C. 
PROGRAM LISTINGS 


A. RESONATOR FREQUENCY 


program SSMSL 
This program calculates resonator frequency 


QO 


integer i, n, m, compare, count 

real*8 a, w, 1, d, t, h, Er, f 

real*§ kn, beta, omega, rl 

real*8 pi, EO, Mu0, c 

real*8 gamls, gam2s 

real*8 5zin, Jzkn, Jzbe, Jz1, Jxkn, Jxbe, Jx1 
real*8 cmp, limit, odlf, newf, limf 


logical llog 


complex*16 gamal, gama2, gama3, Ctl, Ct2, Ct3 
complex*16 Zen, Zenr, Zed, Zer, Zhn, Zhd, Zhnr, Zhr, Zden 
complex*16 Ze, Zh, Zzz, Zzx, Zxz, Zxx 

complex*16 K11, K12, K22, j, skl1, sk12, sk22, det 
character*8 fname 

character*14 sname 





open(2, file=fname 
‘plt' 


open(3, file=sname) 
write.) ' file name 


' fname 
write(3,*) ' file name 


‘ sname 


io] 


c parameter 





(ds 
EO = 1.e-9/(36.*pi) | 
Mu0 = 4.*pi*l.e-7 
limit = 1.e-2 oa 
compare = 1] | 
count = 1 ! 
llog = .true, 1 











22 


OOOO QQ 


RO} 


‘ Jzin = cos(kn*w) — 2. ae UA 





call READF(w,1,a,d,t,h,Er) 


write(2,22) w,],a,d,t,h,Er 


format, "Unit : milimiter' 


,/,t5,'Micro strip line width and length : ',2(e12.6,3x) 


& __ ,/,t5,'Rasonator width : ',e12.6 
& _ ,/,t5, iheights of layers i, 2 and 3 : ',3(e12.6,3x) 
& 145, 'Permittivity of substrate Er 47.3) 


Decide propagation constant Beta 
Applying a root — seeking process 


Calculate the summation 


print*, 'Enter the number of outer summation loop' 


read*, n 
n= 20 
write(2,*) 'The number of outer summation loop = ')n 
write(3,*) 'The number of outer summation loop = ‘jn 


pie ‘Enter the number of inner summation loop' 
read*, 

write (9. * 'The number of inner summation Icop = wm 
write 3 'The number of inner summation loop = 'jm 
rl = 10.41 


initial value 

print”, ‘Enter initial frequency GHz' 
read*, f 

f = f*1.e9 


print*, Frequency Ee | 
write(2, *) 'Frequency = 'f 


omega = 2.*pi*f 





K11 = cmp]x(0.,0. 
K12 = cmplx(0.,0. 
K22 = cmplx(0.,0. 
do 20 i = 1, 


kn = (i - 5)" pi/a 


+ 2.*(1-cos(kn*w Lea he 
Jzkn = 2.*sin katy (knt 3.*Jzin/(kn*w)**2 
Jxkn = 2*pi*sin(kn*w)/(pi**2 - (kn*w)**2) 


Compute summation loop for beta 
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skl1 = cmplx(0.,0. 

sk12 = cmplx(0.,0. 

sk22 = cmplx(0.,0. 
do 40 im = 1 


beta = (real(im) — .5)*pi/rl 

Calculate the gammal and gamma3 

gamls = kn**2 + beta**2 — E0*Mu0*omega**2 
call gamct(gamal,gam1s,Ct1,d) 

call gamct(gama3,gam1s,Ct3,h) 

gam2s = kn**2 + beta**2 — E0*Mu0*Er*omega**2 
call gamct(gama2,gam2s,Ct2,t) 

Calculate Ze 


Zen = gama2*Ct1/Er + gamal*Ct2 
Zenr = -j/(omega*E0) 


Zed = Ctl*Ct2 + Ct1*Ct3*gama2/gama3/Er 
+ Ct2*Ct3 + Er*gamal /gama2 

Ze = Zenr*Zen/Zed 

Calculate Zh 

Zhn = gamal*Ctl + gama2*Ct2 

Zhd = gamal*Ctl*gama2*Ct2 + gamal*Ctl*gama3*Ct3 
+ gama2*Ct2*gama3*Ct3 + gama2**2 

Zhnr = j*omega*Mu0 

Zh = Zhnr*Zhn/Zhd 


Calculate 


Zden = kn**2 + beta**2 

Zzz = —(Ze*beta**2 + Zh*kn**2)/Zden 
Zzx = —kn*beta*(Ze — Zh)/Zden 

Zxz = Z2x 

Zxx = —(Ze*kn**2 + Zh*beta**2)/Zden 


Calculate Jz1 
Jzbe = 4.*pi*cos(beta*l)/(pi**2 — (2.*beta*l)**2) 


Jz) 


= Jzkn*Jzbe 
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Calculate Jx1 


Jxbe = eeipee Die) — sin(beta*l)/(beta*l)**2 


Jxl = Jxkn*IJxbe 
Calculate Matrix coefficient 
skl1 + Jz1*Zzz*Jz1 


skl12 + Jz1*Zzx*Jx1 
sk22 = sk22 + Jx1*Zxx*Jx1 


fa 
~~ 
— 
_ 
eu 


continue 


Kll = Kll + skll 
K12 = K12 + sk12 
K22 = K22 + sk22 


continue 


det = K11*K22 — K12*K12 
print*, 'determinent = ',det 
write(2,*) 'determinent = ',det 


call fdsort(f,det,count) 
count = count + 1 


if(llog) then 


print*, 'Make sure opposite sign between new and old one' 
print*, 'Not opposite sign then Enter 9 ' 
print, a Peo sign then Enter 0! 


itGikeeq),th 
if(ijk.eq.9) then 
oldf = f 


0 to 1 
endif 
log = .false. 
endif 


cmp = real(det) 
limf = tat — oldf) 


print*, 'check limit freq ',limf 
write(2, *) ' check limit freq = 'limf 
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subroutine slval(fun,oxv,xv,newxv,sn) 


integer sn 
real*8 fun,oxv,xv,newxv,pOSxv,negxv 


if(sn_.eq. 1) then 


if(fun .gt. 0) then 


posxv = xv 
negxv = OXxv 
else 
negxv = Xv 
posxv = OXv 
endif 
newxv = (posxv + negxv)/2. 
sn = sn + 2 
else 
if(fun .gt. 0) then 
posxv = xv 
else 
negxv = Xv 
endi 
newxv = (posxv + negxv)/2. 
endif 
return 
end 


He RR HR RIK KR KK I RK AK A EK KK KR OK RK EK KK 2h 2 aK ak 2 kK KK OK 


subroutine fdsort(fqy,sdet,icount) 

integer icount, ic 

real*8 freqy(100),fgy,rtemp, magdet, tfqy 
complex*16 detmt(100),sdet,ctemp, tsdet 


tfay = fay 
tsdet = sdet 


if(icount .eq. )) then 
freqy(1) = tigy 
detmt(1) = tsdet 
elseif(icount .ne. 100) then 














do 10 i = 1, icount-1 


rtemp = pieYy 
ctemp = detmt(i) 


ss ase .gt. tiqy) then 


freqy(i) = tfay 
tiqy = sien 
detmt(i) = tsdet 
tsdet = ctemp 
endif 
continue 


freqy(icount) = tfqy 
‘detmt(icount) = tsdet 
ic = icount 


elseif(icount .eq. 100) then 
write(3,110) 
format(/ /,t5,'Frequency',t20,'Mag of det',t35,'Determinent') 


do 20 j = ljic 
magdet = ca Sr as -+ aimag(detmt(j))**2) 
write(3,210) freqy( yp smagdet ,detmt(j) 
bial +5, f15.2,t20,{15.3,t35,2(e15.8)) 
continue 
endif 


return 
end 


69 














oO 


71 
10 


Calculate Zh 
Zhn = gamal*Ctl + gama2*Ct2 
Zhd = gamal*Ctl*gama2*Ct2 + gamal*Ct1*gama3*Ct3 


& + gama2*Ct2*gama3*Ct3 + gama2**2 


Zhnr = j*omega*Mu0 
Zh = Zhor*Zhn/Zhd 


trmn = ~(B(m)**2*Ze + kn**2*Zh)*Jzkn**2/(B(m)**2 + kn**2) 


sum = trmn + sum 
continue 
mag(m) = sqrt(real(sum)**2 + aimag(sum)**2) 


print71, m,B(m),mag(m),sum 

write(2,71) m,B(m),mag(m),sum 

format(2x,i2,1x,f12.6,1x,f12.6,1x,2(f12.6,1x)) 
continue 


*itmag(n)cmagtn-2) and, magn). )) th 
if(mag(m).lt.mag(m-1) .and. mag(m).It.mag(m+1)) then 
print*, ‘Suspended stripline eend = age 
write oe 
write(2,*) ‘Suspended stripline BETA = ',B(m) 
be = B(m) 


m 
endif 
continue 


Ereff = (be*c/omega)**2 
print*, 'Ereff = ',Ereff 
write(2,*) 'Ereff = ',Ereff 


dth=d+t+h 

call impd(Z0,Ereff,a,dth,t,w) 
print*, 'Z0 = ',Z0 
write(2,*) 'Z0 = ',Z0 
close(1 

close(2 


write(3,95) f, be, Ereff, Z0 
format(2x,e12.6,2x,f8.3,2x,f9.6,2x,f9.5) 


ge to ll 
close(3) 


stop 
end 














Input data file for frequency 


"bezdata.out' 
'f31' 31 
'f32' 32 
'f35' 35 
'f36! 36 

: 'f39' 39 
'f40' 40 
"41! 41 
'f85' 85; 
'f86' 86 
'f93' 93 
'f94' 94 
'f101' 101 
'f102' 102 





Input data file for strip line demensions 


Data 2. 
microstip line width and length 0.60000 0.00000 
Waveguide width :1.27000 
the neight of layers 1, 2 and 3 0.31800 0.12700 0.19000 
the permittivity of substrate :002.220 

» Data 3. 
microstip line width and length :1.50000 0.00000 
Waveguide width :3.56000 

. the height of layers 1, 2 and 3 :0.76300 0.12700 0.89000 
the permittivity of substrate :002.220 


Output data file 


Frequency Beta Ereff Z0 
310000K+11 683.772 1.109129  91.81764 
320000E+1T 705.829 1.109129  91.81764 
850000E+11 772.001 1.109129 91.81764 
360000E+11 794.058 1.109129 91.81764 
390000E+11 860.229 1.109129  91.81764 
A00000E+11 887.601 1.122532 91.26783 
A10000E+11 909.792 1.122532 91.26783 
850000E+11 1953.919 1.204642  §2.52280 
860000E+11 1976.906 1.204642  §82.52280 
930000E+11 2137.817 1.204642  82.52280 
.940000E+11 2160.805 1.204642 82.52280 

. 101000E+12 2321.716 1.204642 82.52280 | 
.102000E+12 2344.703 1.204642 §2.52280 i 

















Data file for input (frequency, beta, Ereff, Z0) 
{32' 32 31 705.829 683.772 1.109129 91.81764 8 
36 35 794.058 772.001 1.109129 91.81764 8 
'f40' 40 39 887.601 860.229 1.109129 91.81764 8 
f86' 86 85 1976.906 1953.919 1.204642 82.5228 6 SS 
94 93 2160.805 2137.817 1.204642 82.5228 6 
'f02' 102 101 2344.703 2321.716 1.204642 82.5228 6 


3 
a 


So 
= 


Data file for input (resonance stub length) 
‘gpl’ 3.280426 3.423938 

'gp2' 3.425125 3.567502 

'gp3' 3.511975 3.655768 

"gp4' 3.567063 3.711716 

'gp5' 3.604022 3.748665 

'gp6' 3.629065 3.773567 

‘gp7' 3.646109 3.790344 

'sp8' 3.657755 3.801590 

'gpl' 2.789856 2.900597 

'gp2!' 2.929428 3.042618 

'gp3' 3.013024 3.128288 

'gp4' 3.067654 3.181847 

'gp5' 3.104645 3.218735 

'gp6' 3.129601 3.244049 

'gp7' 3.146843 3.261261 

tep8! 3.158684 3.272965 

'gp1' 2.398560 2.490021 

'gp2' 2.529320 2.622116 

'gp3' 2.611054 2.704331 

'gp4' 2.664520 2.758667 

'gp5' 2.701723 2.794312 

'gp6' 2.726334 2.819634 

'gp7' 2.743927 2.837242 

'gp8' 2.756226 2.849289 

'g05' 1.211647 1.230606 

'g10' 1.259331 1.278259 

"g15' 1.285843 1.304749 ; 
'g20' 1.300919 1.319977 ‘ 
"g25' 1.310349 1.329392 : 
'g30' 1.316086 1.334999 
'g05' 1.075264 1.090897 ‘| 
'g10' 1.121864 1.137848 | 
'g15' 1.147645 1.163574 \ 
'g20' 1.163193 1.179146 

'g25' 1.172264 1.188179 

'g30' 1.178024 1.193985 

'g05' 0.959969 0.973572 

'g10' 1.005066 1.018852 ; 
'g15' 1.030686 1.044350 . 
'g20' 1.045738 1.059464 
'g25' 1.055008 1.068817 
'g30' 1.060890 1.074570 
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output data file for gap and parallel capacitance 
Freq(GHz) Gap/2 Cgap({F) Cpar(fF) 
.3200E+11 .1000E-03 1.4789 55.8182 
3200E+11 .2000E-03 1.0326 45.8325 
.3200E+11 .3000E-03 8946 40.4891 
.3200E+11 .4000E-03 .8200 37.3474 
.3200E+11 .5000E-03 .7594 35.3625 
.3200E+11 .6000E-03 .7176 34.0614 
.3200E+11 .7000E-03  .6859 33.2005 
.3200E+11 .8000E-03 .6594 32.6301 
.3600E+11 .1000E-03 1.5487 60.9670 
.3600E+11  .2000E-03 1.1390 48.8567 
.3600E+11 .3000E-03 .9906 42.7917 
.3600E+11 .4000E-03 .8371 39.3528 
.3600E+11 .5000E-03 .7638 37.1171 
.3600E+11 .6000E-03 .7307 35.6390 
.3600E+11 .7000E-03 .7012 34.6583 
.3600E+11 .8000E-03 .6784 34.0021 
A000E+11  .1000E-03 9855 67.4308 
4000E+11 .2000E-03 .5812 53.0624 
4000E+11 .3000E-03  .4080 46.0020 
AQ00E+11 .4000E-03 .3403 41.8948 
A000E+11 .5000E-03 .2266 39.3973 
A000E+11 .6000E-03 .2156 37.7064 
AQ00E+11 .7000E-03 .1925 36.5680 
A000E+11 .8000E-03 .1688 35.8057 
S600E+11 .5000E-04 .1896 20.3573 
8600E+11 .1000E-03 .1417 16.8361 
8600E+11 .1500E-03 .1209 15.0875 
8600E+11 .2000E-03 .1154 14.1369 
8600E+11 .2500E-03 .1090 13.5670 
8600E+11 .3000E-03 .1017 13.2336 
.9400E+11  .5000E-04 1893 21.5289 
9400E+11  .1000E-03 1504 17.6024 
9400E+11 .1500E-03 .1257 °15.7299 
.9400E+11  .2000E-03 1148 14.6790 
9400E+11 .2500E-03 .1073 14.08 3 
.9400E+11  .3000E-03 1050 = 13:7098 
.1020E+12 .5000E-04 .2249. 22.9300 
.1020E+12 .1000E-03 1667 18,8290 | 
.1020E+12  .1500E-03 1346 16.4805 
.1020E+12 .2000E-03 1236 15.3471 
.1020E+12 .2500E-03 .1190 14.6777 
.1020E+12 .3000E-03 1101 14.2771 | 
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